Nonlocal Charge Transport Mediated by Spin Diffusion in the Spin-Hall Effect Regime 
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A nonlocal electric response in the spin-Hall regime, resulting from spin diffusion mediating charge 
conduction, is predicted. The spin-mediated transport stands out due to its long-range character, 
and can give dominant contribution to nonlocal resistance. The characteristic range of nonlocality, 
set by the spin diffusion length, can be large enough to allow detection of this effect in materials such 
as GaAs despite its small magnitude. The detection is facilitated by a characteristic nonmonotonic 
dependence of transresistance on the external magnetic field, exhibiting sign changes and decay. 



The spin Hall effect (SHE) is a phenomenon arising 
due to spin-orbit coupling in which charge current pass- 
ing through a sample leads to spin transport in the trans- 
verse direction [l[ . This phenomenon has been attracting 
continuous interest, partially because of the rich physics 
and diversity of SHE mechanisms 0] and partially be- 
cause SHE enables generating spin polarization on a mi- 
cron scale and electrical detection of spin-polarized cur- 
rents, which are the key ingredients of spintronics [3(. 
Theoretically, two main types of SHE have been studied, 
(i) extrinsic, being due to spin-dependent scattering on 
impurities [l| and (ii) intrinsic, arising from the spin-orbit 
terms in the band Hamiltonian 0] . Both extrinsic and in- 
trinsic SHE have been detected experimentally 0, 0, Q 
using optical techniques. Reciprocal SHE (that is, trans- 
verse voltage induced by a spin-polarized current) was 
observed in Al nanowire where ferromagnetic con- 
tacts were used to inject spin-polarized current into the 
sample, and in Pt film [lOj. 

Here we show that the SHE relation between charge 
current and spin current leads to an interesting spin- 
mediated nonlocal charge transport, in which spins gen- 
erated by SHE diffuse through the sample and, by reverse 
SHE, induce electric current elsewhere. The range of non- 
locality of this charge transport mechanism is of the or- 
der of the spin diffusion length £ s = y/D s T s , where D s is 
the spin diffusion constant, and r s is the spin relaxation 
time. The observation of such nonlocal charge transport 
due to SHE can be made fully electrically, which rep- 
resents a distinct advantage compared to the methods 
relying on the sources of spin-polarized current 0, [l(| • 
Although the nonlocal electrical signal, estimated below, 
is small, by optimizing the multiterminal geometry one 
can enhance the nonlocal character of the effect and eas- 
ily distinguish it from the ohmic transport. 

The main distinction of the spin-mediated electrical 
effect considered in the present work from related ideas 
discussed earlier 11, HjIlH 14 1, is that here we identify 
a situation in which the spin-mediated charge transport, 
due to its nonlocality, dominates over the ohmic contri- 
bution. In particular, Refs. 11, 3] considered a correc- 



tion to the bulk conductivity resulting from spin diffusion 



FIG. 1: Nonlocal spin-mediated charge transport schematic. 
Charge current j c applied across a narrow strip generates, via 
SHE, a longitudinal spin current j s - After diffusing over a 
distance x ~ l s 3> w, the spins induce a transverse voltage 
V12 on the probes 1 and 2 via the reciprocal SHE. For a nar- 
row strip, w <^ is, the spin-mediated nonlocal contribution 
exceeds the ohmic contribution that decays as e — w \"\/ w m 



and SHE, which was small in magnitude and thus could 
only be made detectable by its characteristic contribution 
to magnetoresistance. In Ref.[l3| a spin-assisted electri- 
cal response was predicted in a multiterminal mesoscopic 
system, studied numerically in a weak disorder regime 
with the mean free path comparable to the system size. 
In this case the spin-dependent fraction of the multiter- 
minal Landauer-Buttiker conductance is also small. 

The geometry that will be of interest to us is a strip of 
width w narrow compared to the spin diffusion length £ s , 
with current and voltage probes attached as illustrated 
in Fig[TJ The nonlocal charge transport manifests itself 
in a voltage Via across the strip measured at distances 
x ~ l s from the current source (S) and drain (D). Since in 
the purely ohmic, non-SHE regime the voltage away from 
the source decays exponentially on the length scale w/tt, 
our nonlocal spin-dependent voltage can be easily made 
to exceed the ohmic spin-independent contribution. Fur- 
thermore, the spin-mediated effect will exhibit a charac- 
teristic oscillatory dependence on the in-plane magnetic 
field arising due to spin precession during transport. At a 
distance |ac| ~ £ s from the source it will oscillate and de- 
cay as a function of magnetic field on the scale ojb ~ 1 /t s , 
where is the electron spin Larmor frequency. 

We consider, as a simplest case, an infinite narrow strip 

-co < x < +co, —w/2 < y < to/2, (m < i s ), 
as illustrated in FigQ] We assume, without loss of gen- 
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erality, that current source and drain leads are narrow, 
connected to the sample at the points (0,±^w). 

The electric potential in such a sample, described by 
ohmic conductivity <r, can be found as a solution of 
the 2D Laplace's equation with the boundary condition 
j y (x,y = ±|w) — IS(x), where I is the external current. 
Solving it by the Fourier method, we find 



dk 



2nak cosh(ifcw) 



(1) 



In what follows we will need the electric field tangential 
component at the boundaries of the strip, E Xi ±(x) = 
—d x tp±(x), where + and — signs correspond to the strip 
upper edge (y — +w/2) and lower edge (y — —w/2). 
This component of the electric field is found from the 
potential (JTJ) at the boundaries of the sample as 



E x ,±(x) = T 



2/signx 



i\x\ 



oddn>0 



=F 



wa sinh x ' 



(2) 



where x — ttx/w. Potential difference between the strip 
edges, Aip — ip+(x) — (p-(x), evaluated in a similar way, 
decreases as e~^\ x \l w at \x\ > w. 

Below we focus on the case of extrinsic SHE, when the 
fc-linear Dresselhaus and Rashba terms in electron spec- 
trum are negligible. Then the spin fluctuation generated 
by SHE evolves according to the diffusion equation, 

D s d 2 s(x, y) - E(x,y) - — s(x,y) = 0, (3) 

where s is the z-component of the spin density, D s is the 
spin diffusion coefficient, r s is the spin relaxation time. 
The source term S in Eq. ([3]) describes spin current arising 
due to the spin Hall effect: 



-■y(x,y) = V.j s = d a {fl a e a p 1 Ep(x,y)), 



(4) 



where S is the spin Hall conductivity. In the presence 
of the fc-linear Dresselhaus and/or Rashba interaction, 
spin transport is more complicated due to SO-induced 
precession and dephasing [HI, [l(| . The modification of 
the nonlocal electric effect in this case will be briefly dis- 
cussed at the end of the paper. 

Since V x E = 0, the spin current source ([4]) vanishes 
in the bulk and is only non-zero at the strip boundaries: 

E(x,y) = (3 s 8{y-\w)E x , + {x)-{3 s 8{y+\w)E x ,_{x). (5) 

The distinction from the ohmic contribution becomes 
most clear when our strip is relatively narrow, ai<! s . 
In this case the spin current and spin density induced by 
SHE are approximately constant across the strip. Thus 
we can integrate over y and solve a one-dimensional 
spin diffusion problem. Suppressing the y dependence 
in Eq.©, we take E(ir) = (3 s E x . + (x) — j3 s E x -{x). 



Solution of Eq.© in the Fourier representation reads: 
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Pk 
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Dsk2 + 1/Ts , P k = ^l dxZ(x)e^.(6) 



This expression can be simplified by noting that H(x) is 
an odd function of x and that the integral over x con- 
verges at | a; | w, while we are interested in the harmon- 
ics with much lower fc ~ l/£ s *C 1/w. Sending fc to zero 
in the integral, we obtain 



Sk 



1 



ikG 



7r D s k 2 + 1/t s 
where the spin dipole G is given by 
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G = 'E.(x)x dx 
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(7) 



(8) 



(we used Eqs.([2|),((5]) to evaluate this expression). 
Now we can find the spin current 

J s {x) = -D s d x s(x), 



(9) 



where the spin density s(x) is obtained by the inverse 
Fourier transform of Eq.([7]). Using Eq.©, we find 



J s (x) 



2a -L ' 



\x\/i B 



(10) 



(this expression is valid for x not too close to the source, 
| a: | ^ w). The expression (|I0[) gives the net spin current 
rather than the spin current density, as we have been 
solving a ID diffusion problem. This spin current gener- 
ates a voltage across the sample 

6V( X ) = & J 'M = ^fe-W., (11) 
a 2£ s a z 

where f3 c describes charge current arising in response to 
spin current, j" — (3 c e a pj^ . Relating (3 C to the spin 
Hall conductivity as /3 C = Ps/<r, we write the nonlocal 
response (jlll) as a transresistance 



R n \(x) 



SV(x) 1 (13, 



HL e -\*\/i. 

ah 



(12) 



We emphasize that for the extrinsic SHE |l[ , the spin 
current is established on the length scale of the order of 
the electron mean free path £, here taken to be much 
smaller than the strip width w. Thus the inhomogeneity 
of the charge current j c (FigQ]) on the length scale set by 
w does not affect our analysis. The estimate (fT2"j) for the 
nonlocal voltage is therefore accurate as long asi«>l 
We now compare the magnitude of the nonlocal con- 
tribution IT2|) for several materials where extrinsic SHE 
has been observed. For the transresistance IT2|) to be 
large, one would like to have a material with a large ratio 
f3 s /a, and a large spin diffusion length £ s . For Si-doped 
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GaAs with electron density n = 3 x 10 16 cm -3 , the 3D 
charge conductivity, spin Hall conductivity, and spin dif- 
fusion length, reported in 0,0], are given by a^jj w 2.5 x 



71m 



P. 



s3D 



5 x lO- 7 ^- 1 



/im 



9 /im. 



Our two-dimensional quantities a, (3 S are related to the 



3D quantities as a = <J3dw z , (3 s = /3, 



S3D' 



where 



is the sample thickness. Taking w z = 2/im[5|, 0], and 
choosing the sample width to be w = 0.5 /im, we esti- 
mate the transresistance (fT2l) as 



Rnl(x) 



1(T 



x e 



[Ohm] 



(13) 



Although small, it by far exceeds the ohmic conduction 
contribution which at a distance x is proportional to 
a~ 1 e~ n > x >' w . Indeed, for a typical x ~ £ s the ohmic con- 
tribution is negligibly small: e^V™ ~ e -57 K iQ-2i.s 

In the case of InGaAs, the 3D charge conductiv- 
ity and 3D spin Hall conductivity have values similar 
to those quoted above for GaAs (see Ref. Q), (73D ~ 
2.5 x lO^fi-VaT 1 , (3 S3D « 5 x 10~ 7 Q^fiuT 1 , while 
spin diffusion length is considerably shorter, £ s w 2 /im. 
Therefore, in order for the nonlocal voltage (fT2|) at 
|x| ~ t„ to exceed the ohmic contribution, propor- 



tional tO <7 



-tt\x\/w 



the sample width w must satisfy 



w < 7r4/(21n((j 3D //33D)) « 360 nm. 

Another material exhibiting extrinsic SHE is ZnSe Q . 
For carrier concentration n — 9x 10 18 cm -3 the 3D charge 
and spin Hall conductivities are given by 03 d ~ 2 x 
10 _1 fi" 1 Mm- 1 , j3 s3D w 3 x lO^fi -1 ^™ -1 , having the 
ratio /3s3d/c3_d ~ 1.5 x 10 -5 about ten times smaller 
than in GaAs and InAs. The spin diffusion length in this 
material is comparable to that in InAs, l s ~ 2 /im. 

Extrinsic SHE has been also demonstrated in metals, 
Al (see Ref. 0) and Pt (see Ref. (Hill)- In A1 > P»sd ~ 
3 x 10 _3 O _1 J um _1 , the ratio of spin Hall and charge con- 
ductivities is Ps 3 ul G 'iD ~ 1 x 10 -4 , while the spin dif- 
fusion length is £ s w 0.5 /im. Therefore, to separate the 
spin effect from the ohmic contribution, one needs to fab- 
ricate samples with w <C ir£ s /(2\n(<j3D/ 03d)) ~ 85 nm. 
Although the ratio Pssd/vzd ~ 0.37 is large in Pt, 
the observation of the nonlocal effect in this material 
is hindered by its extremely small spin diffusion lcngh, 
l s « 10 nm |19j. We therefore conclude that GaAs sys- 
tems seem to provide an optimal combination of param- 
eter values for the observation of the nonlocal transport. 

We now analyze the effect of an in-plane magnetic field 
on the transresistance (ITBl . In the presence of magnetic 
field, the spin diffusion equation ([3]) is modified as 



D s d 2 s - 3 



h[u B xs] = 0, 



(14) 



where lub = gfJ-B^ is the Larmor precession frequency, 
fiB = 9.27 x 10~ 24 J/T is the Bohr magneton and g is 
the ^-factor. As we shall see below, the interesting field 
range is ^ D s /w 2 . Since in this case the variation 
of spin polarization across the strip is negligible, we can 
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FIG. 2: Nonlocal resistance 7? n i = SV(x)/I sd , Eq.([I5)). in 
units of Ro = (3 c f3 a w /2a 2 l a vs. the in-plane magnetic field 
for several values of x (see Figfl}. [Here lob and t s are the 
electron spin Larmor precession frequency lub and dephasing 
time.] The nonlocal response increases at weak fields, lobt s J; 
£ s /\x\, changes sign at lubt s ~ ^ s /|a;[, and is suppressed at 



UJBTi 



> 



\x\, simultaneously exhibiting oscillations. 



again integrate over the y coordinate and solve a one- 
dimensional diffusion problem. 

For the magnetic field parallel to the x axis, Eq.fl 
takes the following form in the Fourier representation 



g(k) 
g(k) uj b 
-ujb g(k) 

















- 

















(15) 



where g(k) — D s k 2 + l/r s . For the situation of interest, 
when only the z component of the source S is nonzero, 
the solution for s z is given by 



~i(Dsk 2 + l/T s ) 
(D s k 2 + l/T s ) 2 +u 2 



(16) 



Following the same steps as in the absence of magnetic 
field (notice that the source term E z is not affected by 
the magnetic field), we obtain the spin current, 



J s(x) = ^^Re 

2(7 



q+e 



-9+M 



(17) 



where q+ — yl + iujbt s /£ s . 

The nonlocal response due to the voltage induced by 
the spin current (117p . found as SV(x) = ^J s (x), is 



Roi(x) = 



SVjx) = f3 c p s w 
I ~ 2a 2 



Re 



q+e 



-g+m 



(18) 



This expression is simplified in the limit of strong mag- 
netic field, lobt s 3> 1, by factoring an oscillatory term: 



TT 

4 I£ 



(19) 
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where r\ = \Jojbt s /2. Compared to the result found in 
the absence of magnetic field, Eq. fTTj) . the nonlocal volt- 
age SV(x) is amplified by a factor of V2r], decaying on 
a somewhat shorter length scale I = £ s /rj. 

The dependence of R n \, Eq. fTg)) , on the in-plane mag- 
netic field is illustrated in FigJ5] Enhancement of 5V at 
weak fields, wbt s < 1, is followed by a sign change at 
ubt s ~ 1, and suppression at u>bt s > 1. The zeros of SV 
can be found approximately for £ s > \x\ from Eq. (fl~9l) : 



to B , n T s w2^ 2 (n-l/4) 2 4/|x|, 

with integer n > 0. (The condition ^ s /|a;| > 1 ensures 
that io n T s 1, necessary for Eq. flT))) to be valid.) 

For GaAs, g = -0.44 and t s ~ 10ns (see Ref.0]), and 
therefore the field necessary to observe the oscillations 
and suppression of R n \ at ~ £ s , is quite weak: 



In conclusion, spin diffusion in the SHE regime can give 
rise to nonlocal charge conductivity. A relatively large 
nonlocality scale, set by the spin diffusion length, can be 
used to separate the spin-mediated transresistance from 
the ohmic conduction effect. In a narrow strip geometry, 
the transresistance has a nonmonotonic dependence on 
the external in-plane magnetic field, exhibiting multiple 
sign changes and damping. Our estimates indicate that 
observation of the nonlocal conductivity is possible for 
currently available n-doped GaAs samples. 
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and H.-A. Engel. This work is supported by NSF MR- 
SEC (grant DMR 02132802), NSF grants DMR 0541988 
and PHY 0646094, and DOE (contract DEAC 02-98 CH 
10886). 



B ~ B* — 



2mT. 



(20) 



will change 
., decreasing 



5MB t s 

The transresistance measured at \x\ = £ s 
sign at the fields B w 11.1 mT, 60.4 mT, .. 
in magnitude as illustrated in Figj^l 

Nonlocal electric transport can result not only from 
the extrinsic spin scattering mechanisms discussed above 
but also from the intrinsic spin-orbital effects. Below 
we briefly consider this effect and present a rough es- 
timate for a 2D electron gas with Rashba spin-orbit 
coupling, Hso — az.[cr x p], where cr and p are elec- 
tron spin and momentum, z is the unit normal vector, 
and a is SO interaction constant (c/. Ref. 13].) Po- 
tential scattering by impurities leads to Dyakonov-Perel 
spin relaxation with spin diffusion length £ s = h/m*a, 
where m* is the effective electron mass. In such a sys- 
tem, unlike extrinsic SHE, electric field induces spin den- 
sity rather than spin current 11 , [23| ■ The in-plane spin 
polarization s cx z x j c , induced by the source-drain 
current, will diffuse along the strip (FigfTJ), forming a 



profile s(x) ~ am*re" 



'I/a. (For an estimate we 



used the spin diffusion equation derived in [15] , Eq.(13), 
which has a form similar to Eq. (|14[) with a source term 
3(x, y) ~ a 3 (pFT) 2 em*a~ 1 zx j c (x, y), where pf is Fermi 
momentum, and r is the momentum relaxation time.) 

By Onsager relation, the spin density creates electric 
current, j' ~ eorppTZ x s, giving rise to transresistance 



-R n i(x) 



Tit 



w 



(21) 



For a GaAs quantum well with electron density n = 
10 12 cm -2 , mobility fi = 10 5 cm 2 /Vs, spin orbit split- 
ting apF = 100/iV, and width w — 0.5/im, we esti- 
mate R n \ ~ 10 _5 Ohm, which is somewhat larger than 
the ohmic contribution ~ 10 _6 Ohm at x = £ s 3/Ltm. 
However, larger values of the mean free path in quantum 
wells ~ l//m enhance nonlocality of the ohmic contri- 
bution. This may hinder observation of spin-mediated 
transport despite a larger value of R n \ ( cf. Ref. [l]| ) . 



[1 
[2. 

[3; 
[5 

ir 

[§ 

[9 
[10 

[11 

[12 
[13 

[14 
[15 

[16 
[17 

[is: 

[19 

[20 



M. I. D'yakonov and V. I. Perel, JETP Lett. 13, 467 
(1971). 

H.-A. Engel, E. I. Rashba, and B. I. Halperin, in: Hand- 
book of Magnetism and Advanced Magnetic Materials, 
vol. V (Wiley, 2006). 

S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. 
M. Daughton, S. von Molnar, M. L. Roukes, A. Y. 
Chtchelkanova, and D. M. Treger, Science 294, 1488 
(2001). 

J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jung- 
wirth, and A. H. MacDonald, Phys. Rev. Lett. 92, 126603 
(2004). 

Y. K. Kato, R. C. Meyers, A. C. Gossard, D. D. 
Awschalom, Science 306, 1910 (2004). 
J. Wunderlich, B. Kaestner, J. Sinova, and T. Jungwirth, 
Phys. Rev. Lett. 94, 047204 (2005). 

V. Sih, W. H. Lau, R. C. Meyers, V. R. Horowitz, A. 
C. Gossard, and D. D. Awschalom, Phys. Rev. Lett. 97, 
096605 (2006). 

N. P. Stern, S. Ghosh, G. Xiang, M. Zhu, N. Samarth, 
and D. D. Awschalom, Phys. Rev. Lett. 97, 126603 
(2006). 

S. O. Valenzuela and M. Tinkham, Nature 442, 176 
(2006). 

E. Saitoh, M. Ueda, H. Miyajima, and G. Tatara, Appl. 
Phys. Lett. 88, 182509 (2006). 

L. S. Levitov, Y. V. Nazarov, and G. M. Eliashberg, Sov. 

Phys. JETP 61, 133 (1985). 

J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999). 

E. M. Hankiewicz, L. W. Molenkamp, T. Jungwirth, and 

J. Sinova, Phys. Rev. B 70, 24130 1(R) (2004). 

M. I. Dyakonov. larXiv:0705. 27381 

E. G. Mishchenko, A. V. Shytov, and B. I. Halperin, 
Phys. Rev. Lett. 93, 226602 (2004). 

W.-K. Tse and S. Das Sarma, Phys. Rev. B 74, 245309 
(2006). 

T. Kimura, Y. Otani, T. Sato, S. Takahashi, and S. 
Maekawa, Phys. Rev. Lett. 98, 156601 (2007). 

G. Y. Guo, S. Murakami, T.-W. Chen, and N. Nagaosa, 
larXiv:0705.0409l (unpublished). 

H. Kurt, R. Loloee, K. Eid, W. P. Pratt, and J. Bass, 
Appl. Phys. Lett. 81, 4787 (2002). 

V. M. Edelstein, Solid State Comm. 73, 233 (1990). 



